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Abstract. We investigate the dynamics of a large anisotropic spin whose easy- 
axis component is coupled to a bosonic bath with a spectral function J{lu) oc lu". 
Such a spin complex might be realized in a single-molecular magnet. Using the 
non-perturbative renormalization group, we calculate the line of quantum-phase 
transitions in the sub-ohmic regime (s < 1). These quantum-phase transitions 
only occur for integer spin J. For half-integer J, the low temperature fixed-point 
is identical to the fixed-point of the spin-boson model without quantum-tunneling 
between the two levels. Short-time coherent oscillations in the spin decay prevail 
even into the localized phase in the sub-ohmic regime. The influence of the 
reorganization energy and the recurrence time on the decoherence in the absence 
of quantum-tunneling is discussed. 
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1. Introduction 

Understanding the influence of the environment on the non-equilibrium dynamics of 
quantum systems remains one of the challenging questions of theoretical physics. A 
finite number of quantum mechanical degrees of freedom, a large spin or a qubit 
interacting with a infinitely large bath of non-interacting bosons with a continuous 
energy spectrum represents a typical class of model examples for such systems. Its 
simplest version, the spin-boson model [1], has contributed tremendously to our 
understanding of dissipation in quantum systems [2] . 

Single-molecular magnets (SMM) based on Mni2-actetate or Feg-complexes 
behave as large single spins at low temperatures. The magnetic anisotropy of the 
molecules prevent a simple relaxation of the spin. The resulting hysteretic behavior in 
the magnetization has been subject to intensive experimental [3l|4l|5l[6] and theoretical 
investigations [II [3 [9l [1^ (see Ref. for an overview). Using Wilson's numerical 
renormalization group (NRG) method [TTl [T^ , it was shown [T^ that the combination 
of quantum tunneling and antiferromagnetic coupling to a metallic substrate induces 
a Kondo effect in such SMM. The recently developed time-dependent numerical 
renormalization group method (TD-NRG) [TS] has been employed to analyze the 
real-time dynamics of such large spin [16j . 

In this paper we will investigate the equilibrium and non-equilibrium spin 
dynamics of a anisotropic spin coupled to a dissipative bosonic bath. Such a 
spin complex might be realized in a single-molecular magnet in which the physical 
properties depend on the type of coupling to the environment. Previously, detailed 
expansion of the spin-lattice relaxation [17] in powers of the local spin operator have 
been considered to estimate the spin-relaxation rates [8] by Fermi's golden-rule 
arguments. Vorrath and Brandes have reduced the spin-bath coupling to a linear 
term proportional to the Sz operator [lOj. A Weisskopf-Wigner coupling to the 
environment would describe the interaction with a fluctuating quantized magnetic 
field. This might be relevant for recent dephasing experiments ^ in SMM. Here, we 
will restrict ourselves to diagonally coupled spin-boson model interaction [H [10] . 
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The real-time dynamics in the spin-boson model is usually investigated using the 
noninteracting-bilp approximation [U [T5] or by Bloch-Redfield type of approaches 
[T^ . Recently, it has been demonstrated that the non-perturbative numerical 
renormalization group [1111201 [^Tlll2j and its extension to non-equilibrium dynamics, 
the TD-NRG [Ml [15] , is particularly suitable to access the quantum-critical region of 
the spin-boson model in a sub-ohmic environment [201 EI] and is able to describe the 
real-time dynamics on short-time scales as well as exponentially long-time scales |22j . 

In order to gain a better understanding of the non-equilibrium dynamics, we 
also present the equilibrium properties of the model. We will map out the quantum- 
critical line separating a delocalized phase in weak coupling and a local phase at strong 
coupling to the dissipative sub-ohmic bath. The continuous quantum-phase transitions 
found for < s < 1 and integer spin J resemble the phase diagram previously reported 
for the spin-boson model by Bulla et al. [2^. No phase transition was found for half 
integer spin J. Information on the equilibrium dynamics is provided by the equilibrium 
spin correlation function C{uj). The non-equilibrium spin dynamics of a anisotropic 
spin complex coupled to a dissipative bosonic bath is investigated in response to a 
sudden change of magnetic field close to quantum phase transitions. 



2. Modelling of a large spin in a dissipative environment 

The local Hamiltonian modeling single-molecular magnets [5] 

H — Hloc + Ti-I + Hbath (1) 

consists of a single large spin of size J subject to an easy-axis anisotropy energy A and 
quantum-tunneling terms [S1[T3| B2n, which induce transitions between the eigenstates 
of 5.: 

Hioc = - ASl + B2n (5^" + 5f ) . (2) 

n=l,2 

These quantum-tunneling terms stem from a reduction of the U{\) easy-axis spin 
symmetry to the discrete symmetry group of the SMM molecule. Assuming A \B2n\ 
and a positive A , the energies of states | J^) are located on an inverted energy parabola. 
As depicted in Fig. [1] there is a fundamental difference between integer and half- 
integer values of the spin [5] [13] : the two ground states | ± J) are connected via the 
quantum-tunneling terms i?2n only for integer values of J, while for half-integer values 
the ground states remain disconnected. Without further transition mechanism, a half- 
integer anisotropic spin cannot relax into thermodynamic ground state after switching 
off the spin-polarizing external magnetic field. 

We use a generalization lOJ of the spin-boson model (SBM) 

Hi =S.J2^Abl + b-,) (3) 

9 

Hbath = ^ i^qblbq (4) 
9 

to describe the coupling to a dissipative bath. Here, 6j creates a boson in the mode q 
with energy ujq. The coupling between spin and bosonic bath is completely specified 
by the bath spectral function [1] 

J{cj)^7rJ2>^lS{cj-cjq) . (5) 

9 
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(a) (b) 

Figure 1. Level scheme for (a) spin J = 7/2 and (b) J = 4. The two states 
I ± J) are only connected by the quantum tunneling terms B2 and for integer 
spin J. 



The asymptotic low-temperature behavior is determined by the low-energy part of the 
spectrum. Discarding high-energy details, the standard parametrization is 

J (uj) ^ 27r a Lul~'' w\ < Lu < LJc , s > -I (6) 

where the dimensionless parameter a characterizes the dissipation strength. tOc is a 
cutoff energy and is used as the energy scale throughout the paper. The value s = 1 
corresponds to the case of ohmic dissipation. In the usual spin-boson model (SBM), 
the local spin Hamiltonian 

ttSBM ^ ^ r'7\ 

is parametrized by an energy splitting e and a quantum-tunneling term A where ai 
are Pauli matrices [I]. 

3. Numerical renormalization group 

3.1. Equilibrium NRG for a bosonic environment 

In all our calculations we have used Wilson's numerical renormalization group (NRG) 
method. Wilson's NRG method is a very powerful tool for accurately calculating 
equilibrium properties of quantum impurity models (12] . Originally the NRG was 
invented by Wilson for a fermionic bath to solve the Kondo problem [TTl [23] . The 
method was recently extended to treat quantum impurities coupled to a bosonic bath 
pUl [5T] , combination of fermionic and bosonic baths [511 and real-time dynamics out 
of equilibrium [5S1 [Til [TS] . The non-perturbative NRG approach has been successfully 
applied to arbitrary electron-bath couplings strength [501 [211 [HI [22 • 

At the heart of this approach is a logarithmic discretization of the continuous bath, 
controlled by the discretization parameter A > 1; the continuum limit is recovered 
for A — > 1. Using an appropriate unitary transformation [11] . the Hamiltonian is 
mapped onto a semi-infinite tight-binding chain, defined by a sequence of finite-size 
Hamiltonians Hm with the impurity coupled to the open end. The tight-binding 
parameters tm linking consecutive sites of the chain m and m+1 falls off exponentially 
as tm ^ A""*. Each bosonic chain link is viewed as representative of an energy shell 
since its energy Wm also decreases as Wm A~™ establishing an energy hierarchy. 
Both ensure that mode coupling can only occur between neighboring energy shells. 
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which is essential for the application of the renormalization group procedure. To this 
end, the renormalization group transformation R[H] reads 

'Hm+l = R[Hm] 

= AHm + tm (^al^am+l + aL+l^^™) + li'mflL+l^m+l , (8) 

where Tim is the Hamiltonian of a finite chain up to the site m. The annihilation 
(creation) operators of site m are denoted by am {o-ln) "^m represents the energy 
of the bosonic mode of site m. Note that the rescaling of the Hamiltonian Hm by A 
ensures the invariance of the energy spectrum of fixed point Hamiltonians under the 
RG transformation RlTCm]- For a detailed review on this method see Ref. [12\ . 

The RG transformation ([8|) is used to set up and diagonalize iteratively the 
sequence of Hamiltonians Hm- In the first step, only the large spin coupling to the 
single bosonic site to = is considered. It turns out to be sufficient [501 IHl [H] to 
include only the Ni, lowest lying bosonic states, where Ni, takes typical values of 8 — 12. 
The reason for that is a quite subtle one: the coupling between different chain links 
decays exponentially and is restricted to nearest-neighbor coupling by construction, 
both essential for the RG procedure. In each successive step (i) a finite number of 
Nb bosonic states of the next site m + 1 are added, (ii) the Hamiltonian matrices 
are diagonalized and (iii) only the lowest Ns states are retained in each iteration. 
The discarding of high energy states is justified by the Boltzmannian form of the 
equilibrium density operator when the temperature is lowered simultaneously in each 
iteration step to the order Tm oc A~™Wc. 

Denoting the set of low- lying eigenstates by \r)j^ and the corresponding 
eigenvalues Er{N) oc 0(1) at iteration N, the equilibrium density matrix pq is given 

m by 

P0 = ^j:e-^^^\r)^^{r\, (9) 
Zn ^ 

where — J2r^~^^'^ ^^'^ P ^'"^ the order 0(1), such that Tjv = WcA~^/p. 
The thermodynamic expectation value of each local observable O is accessible at each 
temperature Tjv by the trace 

(d)eq = TrUol ^ ^y^e-^''" N{r\d\r)N. (10) 
L J Zn ^ 

The procedure described above turns out to be very accurate because the coupling tm 
between the bosonic sites along the chain are falling off exponentially so that the rest 
of the semi-infinite chain contributes only perturbatively jlll 112] at each iteration to, 
while contributions from the discarded high-energy states are exponentially suppressed 
by the Boltzmann factor. 



3.2. Time- dependent NRG 

While the equilibrium properties are fully determined by the energy spectrum of 
the Hamiltonian, the non-equilibrium dynamics requires two conditions: the initial 
condition encoded in the many-body density operator po and the Hamiltonian Ti-^ 
which governs its time-evolution. For a time-independent Hamiltonian, the density 
operator evolves according to p(t > 0) = e~*^ *Poe*^ *■ AH time-dependent 
expectation values (O) {t) are given by 

(11) 



{d){t) - Tr p{t)0 



= Tr 



m't 



O 
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where set h = I. 

We obtain the density operator po from an independent NRG run using a suitable 
initial Hamiltonian W. For instance, by choosing a large local magnetic field in W, 
we can prepare the system such that Sz is fully polarized. To investigate decoherence 
we choose H}^^ such that 

l^)-^EH (12) 

m 

is an eigenstate of an appropriate Ti-j^^ with the lowest eigenenergy. 

In general, the initial density operator po contains states which are most likely 
superpositions of excited states of . For the calculation of the real-time dynamics of 
the large spin it is therefore not sufficient to take into account only the retained states 
of the Hamiltonian obtained from an NRG procedure. The recently developed 
time-dependent NRG (TD-NRG) [111 [15] circumvents this problem by including 
contributions from all states. It turns out that the set of all discarded states eliminated 
during the NRG procedure form a complete basis set [131 [KJ of the Wilson chain which 
is also an approximate eigenbasis of the Hamiltonian. Using this complete basis, it 
was shown fl4J [15] that Eq. pT|) transforms into the central equation of the TD-NRG 
for the temperature T/y 

N trun 

(0)it) = E E e^'^"-^"'*0™ p-^(m) , (13) 

m— r,s 

where O™^ — (r; to|0|s; to) are the matrix elements of any operator O of the electronic 
subsystem at iteration to, and E",™, E™ are the eigenenergies of the eigenstates |r; to) 
and |s;to) of At each iteration to, the chain is formally partitioned into a 

"system" part on which the Hamiltonian Tim acts exclusively and an environment 
part formed by the bosonic sites to + 1 to A^. Tracing out these environmental degrees 
of freedom e yields the reduced density matrix[T4l [15] 

e 

at iteration m, where po is given by ([9]) using W. The restricted sum X)'' -^Q- 
implies that at least one of the states r and s is discarded at iteration to. Excitations 
involving only kept states contribute at a later iteration and must be excluded from 
the sum. 

As a consequence, all energy shells to contribute to the time evolution: the short 
time dynamics is governed by the high energy states while the long time behavior is 
determined by the low lying excitations. Dephasing and dissipation is encoded in the 
phase factors e*^^""-^™^* as well as the reduced density matrix pl°^{m). 

Discretization of the bath continuum will lead to finite-size oscillations of the real- 
time dynamics around the continuum solution and deviations of expectation values 
from the true equilibrium at long time scales. In order to separate the unphysical 
finite-size oscillations from the true continuum behavior, we average over different 
bath discretization schemes using Oliveira's z-averaging (for details see Refs. [Ml [IS]). 
We average over Nz = 16 different bath discretizations in our calculation. 

Previously, the TD-NRG has been successfully applied to the simple spin-boson 
model [151 HI] , the SMM coupled to a fermionic bath [16] and electron-transfer in a 
dissipative environment [27] . 
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Figure 2. Zero temperature phase diagram in the sub-ohmic regime for J=3. 
The critical coupling strength ac(2J)^ is plotted vs the exponent s in J{u}) for 
A = -0.1, B2 = B4 = 10-3, H = 0. NRG Parameters: A = 2, TVs = 100, Nt = 8. 



4. Equilibrium and quantum-phase transitions 

Integer Spin J and quantum-phase transitions 

In Fig. [51 we plot the zero-temperature phase diagram in the sub-ohmic regime for a 
moderate value of an integer spin (J = 3.) Depicted is the rescaled critical coupling 
strength ac(2J)^ as a function of the bath power-law exponent in the sub-ohmic 
regime. The line marks the quantum-phase transition (QPT) between a localized and 
a delocalized phase. We have obtained the critical coupling strength by investigating 
the NRG level flow which shows that the QPT is associated with a distincted fixed 
point. For s > 1, no QPT exist [201I2T]. 

Analyzing the NRG fixed-point spectra [H] , the same qualitative picture emerges 
for all integer values of J. We noted that the critical coupling strength depends (i) on 
the value of the spin J and (ii) the values of the quantum tunneling rates B2n- Only 
for s ^ 1, the critical coupling becomes independent of and always approaches 

lima,(s) « 1/(2J)2 , (15) 

which can be understood from the analytical form of Tii. 

The fixed-point spectra always agrees with those of the spin-boson model. This 
establishes a mapping of our model at low temperatures onto an effective spin-boson 
model. This mapping becomes intuitively clear by inspecting Fig. [TJb). It shows 
a coupling between the two lowest states | ± J) on the parabola. For T ^ 0, it 
is expected that the model can be replaced by an effective spin-boson model with an 
tunneling rate A which is a complicated function of parameters A and i?2n- Therefore, 
it is not surprising that our calculations for the anisotropic large spin model reveal 
the similarities to the results by Bulla et al. '^T\, who established such a quantum 
phase transition in the SBM. Since the coupling of the two states | ± J) to the bosonic 
bath is proportional to Jz, we find that ac{s){2J)'^ resembles the phase boundary of 
spin-boson model. The critical coupling ac is proportional to 

Such a behaviour is also observed in the equilibrium spin-spin correlation function 
C{uj) = 3m ^ Sz\Sz ^ {u! — iS). The correlation function C{u!) is calculated using 
the sum-rule conserving algorithm |^. We plot C{uj) on a log-log scale for a 
series of coupling constants and the sub-ohmic regime s = 0.5 in Fig. [Sja). The 
high energy feature stemming from the easy-axis energy splitting proportional to A is 
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5 10 15 20 25 30 



Figure 4. NRG level flow of the lowest eigenenergies i?™ for J = 7/2, a = 0.0012 
in the sub-ohmic regime s = 0.5 in comparison with the level flow of the spin-boson 
model for A = and a = 0.05. Parameters as in Fig. |2] 

clearly visible. At lowest coupling strength a = 10"^, the second peak at a; « 3 x 10"'^ 
is well pronounced, originating from the weak effective quantum tunneling rate which 
is a function of parameters of TLioc- A,B2n as well as J. For a < ac, C'(w) vanishes 
according to C{uj) oc uj'^, while C{uj) diverges as w^"* for cj ^ and a > ac- 

4-2. Half-integer Spin J 

A completely different picture emerges for half-integer spins J. Independent of the 
bath exponent s, no phase transition is found. As seen already in Fig. [1] the quantum 
tunneling matrix elements B2 and B4 cannot connect the two states | ± J). The 
model can be mapped onto a spin-boson model with a vanishing tunneling rate A 
at low temperatures. To illustrate this, we compare the NRG level flow |:12 of the 
large spin model with the NRG level flow of a simple spin-boson model and A = 0. 
The fixed-point spectra for m — > 00 is identical for both models as depicted in Fig. [H 
The coupling constant a determines the crossover from high-temperature to the low 
temperature fixed point: this crossover scale accidentally coincides in both models for 




our choice of coupling constants. C{uj) depicted in Fig. ^h) is similar to the results 
of Fig. [Sja) at high energies: the peak due to the easy-axis splitting A is also present. 
However, a pronounced peak at the effective quantum-tunnel energy is absent for any 
half-integer J. 

5. Non-equilibrium dynamics 

For the investigation of the non-equilibrium dynamics, the local spin is prepared in a 
maximal polarized initial state with the expectation value 52(0) = J by applying a 
strong local magnetic field in z-direction. The external magnetic field is switched off 
at i = 0. 

Fig. O displays the real-time dynamics of Sz{t) for a series of coupling strength 
a and an ohmic bath {s = 1), J — 4 and the easy-axis splitting A — 0.1. A weak 
quantum-tunneling was enabled by B2 = 10^^ and B4 ^ 2x 10^'^. As a reference, the 
dynamics of a decouple spin (a — 0) is included in the graph. For small a, damped 
coherent oscillations with a superimposed short-time dynamics of small amplitude are 
clearly visible. The oscillation frequency decreases with increasing a and vanished 
at a finite coupling a. Above that value - here a « 0.075 - only a spin decay is 
observed. This corresponds to the regime 0.5 < a < 1 of the SBM. When approaching 
the critical coupling adJ = 4) « 0.016, Sz{t) only decays insignificantly on a very 
short time scale of the order 0(1) and remains constant for t —>■ 00, as seen in panel 
Fig.[5fa). This behaviour occurs in the localized phase of the model where the effective 
tunneling rate is renormalized to zero. 

A different picture emerges in the sub-ohmic regime. The real-time dynamics 
of a large anisotropic spin with J = 5 is depicted for different coupling strength in 
Fig. [6l The upper panel shows the dynamics on a linear time axis. In the weak 
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Figure 6. Real-time dynamics of tlie spin Sz initially prepared in Sz{0) = J = 5 
for different values of a in tlie sub-ohmic regime T = 10~®. Tlie upper panel (a) 
shows the short time dynamics on a linear time scale, while (b) resolves the real- 
time dynamics on a logarithmic time scale. For couplings a > ~= 0.00061, 
the finite value oi Sz{t ^ oo) is a consequence of the quantum-phase transition 
into a localized phase. Parameters: A = 0.1, B2 = 10"^, B4 = 2 X 10~^, A = 2, 
Ns = 160, Nt = 8, N^ = 16. 

coupling regime, coherent high and low frequency oscillations are clearly visible, in 
addition to a slow decay of the amplitude. At a critical coupling of a = 0.00061, we 
found a quantum-phase transition in the equilibrium NRG calculations. As discussed 
in Sec. 14.11 the quantum-critical point separates a delocalized from a localized phase. 
As seen in Fig. [51 coherent oscillations prevail in the short-time dynamics even in the 
localized phase close to the quantum-phase transition. A finite time is required before 
the strong correlations can build up. In the long-time limit, Sz{t) approaches a finite 
value in the localized phase which should be proportional to the order parameter 
for the QPT. These results are very similar to the dynamics reported in the sub- 
ohmic spin-boson model [52]. In our case, larger number of eigenstatcs of Hioc yields 
additional coherent oscillations absent in the simple model. 

We expect that coherent oscillations prevail to longer time-scales with decreasing 
exponent s. Simultaneously, the critical coupling ac reduces strongly with deceasing 
exponent s, as shown in Fig.[2l 




10"^ 10"' lo" lo' 10" lo' 

t*a 



Figure 7. Sxit) vs t*a for three different coupling strength a = 10~*, 10~^, 10~^ 
for J = 5 (a) and J = 9/2 (b) values of the spin in the s = 0.5 sub-ohmic 
regime. Sx{t) measures the decoherence of the off-diagonal matrix elements of 
pT-ed Parameters: T = 10~^, A = B2n = 0. (c) shows Sx{t) vs t*Q for the same 
values of a as (b) but with A = 2aLJc/s, B2n = 0. NRG parameters as in Fig. [5] 



6. Decoherence 



Decoherence is caused by the interaction of a subsystem with its environment. The 
concept of decoherence [30l [31] provides inside into how a coherence quantum state 
\s) of the subsystem is destroyed by the entanglement with the environment. The 
reduced density matrix [32] 

e 

evolves from a matrix describing such a pure state |s) into a matrix describing an 
ensemble. In an appropriate local basis the off-diagonal matrix elements of p\'^^ 
must vanish for i ^ oo [30l [31] . In a diagonally coupled spin-boson model as described 
by the interaction TLi, Eq. (j4|), this would be the eigenstates of 5*2 ^Oj. In the absence 
of quantum tunneling, Sz also commutes with Tiioc and no energy is exchanged by 
the coupling Tij to the environment. Decoherence is induced by the dephasing with 
continuum of the bath modes each contributing a different phase shift. For J — 1/2, 
this is a well studied problem, and it was shown analytically [33l [34] that the off- 
diagonal component of the reduced density matrix for the spin can be written as 
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Phli^) ~ ^ '"''*VT,i(0)' where r(t) is given by the exact analytic expression 

xn/.x 1 1"°° , T/ N , / \ 1 - cos(uji) ,,„, 
m--l dc.JMcoth(-) J r (17) 

Here T is the temperature. Note that decoherence of the pure quantum state \s) does 
not follow a simple exponential form. Approximating the dephasing by a function 
oc exp{—t/Tcieph), which defines a single "dephasing" time-scale Tdeph, might be 
insufhcient even in very simple models such as the SBM, since T(t) cannot be replaced 
by t/Tdeph- 

We initially prepared the spin system in a pure state \s) comprising of a linear 
combination of all Sz eigenstates 

1^) - 27TT ^ l-> (18) 

m— — J 

by an appropriate choice of the initial Hamiltonian Ti.}^^. 

We can gain some information on the decoherence of a large anisotropic spin by 
measuring the time-dependent expectation value Sx{t) — {Sx){t), which only depend 
on the off-diagonal matrix elements of plj^{t). The quantum-tunneling parameters 
were set to B2n = in the calculations. Fig. [7] shows the dynamics of Sx{t) on a 
rescaled dimensionless time axis ta for a integer (J — 5) and a half-integer (J — 9/2) 
value of the spin in the s = 0.5 sub-ohmic regime. In both cases, we observe an 
oscillatory decay of Sx different from the result for the spin-boson model, Eq. (fT7|) . 
Here, the finite coupling a introduces an effective A via the reorganization energy Eig 

Ei„ = y^ = duj = , s > . (19) 

^ UJq Jo TTtJ S 

The effective energy Em oc {A~ Eig)m'^ of eigenstates |m) of 5*2 induces quantum- 
oscillations whose oscillation frequency is proportional to a for ^ = 0. This a 
dependence of the frequency is revealed by plotting Sx{t) versus the dimensionless 
timescale ta so that the oscillation beats coincide and occur on a time scale of 
ta « 0(1). Since all eigenenergies Em are related to each other by integer values 
proportional to ^ — Eig, the local dynamics is characterized by a recurrence time 
proportional to 1/{A — Eig) at which the original local state would reappear in the 
absence of a spin-environment coupling. However, the envelope function describing the 
time-dependence of the recurrence amplitude decreases on a time scale proportional 
to 1/a. 

We can prolong this recurrence time by an appropriate counter-term in the 
Hamiltonian. The curves of Fig. W[c) are obtained for the same parameters as in 
Fig.[7I^a) and by setting A = laujcj s. Due to the discretization of the bath continuum, 
the cancellation is not perfect, but the oscillation frequency is lowered by almost two 
decades. 



7. Conclusion 



The equilibrium and non-equilibrium dynamics of an anisotropic large spin coupled 
to a dissipative sub-ohmic bosonic environment has been investigated using the non- 
perturbative NRG. Such a large spin might be realized in a single-molecular magnet. 
The coupling of the bath modes to the easy-axis component of the spin enhances the 
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tendency for localization. A quantum-phase transition exists only for integer spin 
values in the sub-ohmic regime, and the quantum-critical point is associated with a 
new fixed point. These results generalized the work of Bulla et al. PU] for the spin- 
boson model. The spin-spin correlation function C{uj) diverges approximately as \uj\~''' 
at and above the critical coupling ac- The critical coupling for s — > 1 is related to the 
critical coupling of the spin-boson model by ac = af^*^/(2J)^. We have shown that 
the fixed point of the half-integer model is identical to the S = 1/2 spin-boson model 
without quantum-tunneling A. No quantum-phase transition is found for half- integer 
J due to the lack of quantum-tunneling between the two states | ± J). 

We investigated the spin decay by switching off an external easy-axis magnetic 
field at t = 0. In the ohmic regime (s = 1), the real-time dynamics of Sz{t) is governed 
by several local frequencies stemming from the easy-axis splitting A and the quantum- 
tunneling matrix elements i?2n- With increasing coupling a, the quantum-tunneling 
is reduced, but the splitting A is increased by the reorganization energy Eig. The 
low-frequency coherent oscillations vanish at some value a but the high-frequency 
coherent oscillations prevail with a very small amplitude. Further increasing of a 
beyond the critical coupling ac drives the system in the localized phase. Here, the 
spin decays weakly on a very short time-scale of the order of the reciprocal cutoff 
I/lOc- a completely different picture emerges in a sub-ohmic environment. Oscillatory 
response is found at short and intermediate time-scales even in the localized phase at 
coupling strength a > ac- 

In the calculation of the decoherence through dephasing, we identified the 
occurring long-time oscillation frequency with recurrence frequency which is dependent 
on the reorganization energy. These oscillations can be reduced and suppressed by 
adding a counter-term in the original Hamiltonian. We have demonstrated that the 
non-perturbative renormalization group approach can describe such equilibrium and 
non-equilibrium dynamics close to a quantum-phase transition. 
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